Classical and quantum pumping in closed systems 
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Pumping of charge (Q) in a closed ring geometry is not quantized even in the strict adiabatic 
limit. The deviation form exact quantization can be related to the Thouless conductance. We use 
the Kubo formalism as a starting point for the calculation of both the dissipative and the adiabatic 
contributions to Q. As an application we bring examples for classical dissipative pumping, classical 
adiabatic pumping, and in particular we make an explicit calculation for quantum pumping in case 
of the simplest pumping device, which is a 3 site lattice model. We make a connection with the 
popular S matrix formalism which has been used to calculate pumping in open systems. 



Pumping of charge in mesoscopic Q and molecular size 
devices is regarded as a major issue in the realization of 
future " quantum circuits" or " quantum gates" , possibly 
for the purpose of " quantum computing" . Of particu- 
lar interest is the possibility to realize a pumping cy- 
cle that transfers exactly one unit of charge per cycle 
[J IE IE 0- I n open systems this "quantization" holds 
only approximately. But it has been argued |(J that the 
deviation from quantization is due to "dissipative" effect, 
and that exact quantization would hold in the strict adi- 
abatic limit, if the system were closed. In this Letter 
we would like to show that the correct picture is quite 
different. In particular we would like to make a proper 
distinction between " dissipative" and " adiabatic" contri- 
butions to the pumping, and to calculate the deviation 
from exact quantization in the latter case. As a starting 
point we adopt the traditional Kubo formula lal , bu t we 
also point out the relation to the " adiabatic" ||J, Ll3 an( i 
to the "S'-matrix" formulations. The present formu- 
lation of the pumping problem has few advantages: It 
is not restricted to the adiabatic regime; It give a "level 
by level" understanding of the pumping process; It al- 
lows the consideration of any type of occupation (not 
necessarily Fermi occupation); It allows future incorpo- 
ration of external environmental influences such as that 
of noise; It regards the "voltage" over the pump as "elec- 
tro motive force" , rather than adopting the conceptually 
complicated view [l2| of having a "chemical potential dif- 
ference". Finally, on the practical level, we give a solu- 
tion for the pumping in a 3 site lattice model. This is 
definitely the simplest pump circuit possible, and we be- 
lieve that it can be realized as a molecular size device. It 
also can be regarded as an approximation for the closed 
geometry version of the two delta potential pump Q • 

The structure of this Letter is as follows: We show 
how to get from the Kubo formalism an expression for 
the pumped charge Q, and explain the distinction be- 
tween "dissipative" and the "adiabatic" contributions. 
Then we give illuminating examples for classical dissipa- 
tive pumping and for classical adiabatic pumping. Next 
we discuss the case of quantum pumping, where the cy- 
cle is around a chain of degeneracies. We show that this 
can be understood as a special case of " adiabatic trans- 
fer" scheme. In order to get a quantitative estimate for 
the pumped charge we consider a 3 site lattice model. 



We get expressions for Q, and express them in terms of 
the Thouless conductance. We conclude by a short dis- 
cussion of the relation between the Kubo formalism, the 
adiabatic formalism, and the S'-matrix formalism. 

Consider a system that has a ring geometry (Fig. la). 
The Hamiltonian is Tl(xi(t), X2(t), x^it)), where x\ and 
X2 are parameters that control the shape of the ring, 
or the height of some barriers, while x^=^=T%(j) is the 
magnetic flux. We use units such that the elementary 
charge is unity. The " generalized forces" are convention- 
ally defined as F k = —dTt/dxk- In particular (F 3 ) is the 
current / through the ring (see remark |!8j1. Consider 
for a moment the time independent Hamiltonian TL{x), 
with x = const, and assume that the system is prepared 
in a stationary state (either pure or mixed). The expec- 
tation value (F h ) of a generalized force is known as the 
"conservative force" or (in case of fc=3) as the "persistent 
current" . The " fluctuations" of the generalized forces are 
conventionally characterized by the real functions: 



C y '(r) 
K 11 {t) 



(i(F l (r)^(0)+^'(0)F l (r))> (1) 



[F*(r),^'(0)]) 



(2) 



Note that both functions have a well defined classical 
limit. Their Fourier transform will be denoted by C tJ (cu) 
and K % HlS) respectively. 

Our interest in the following is in a driving cycle, where 
x = x(t) forms a loop in the 3 dimensional parameter 
space. In linear response theory 8] the non-conservative 
contribution to (-F fc ) is related to x(t) by a causal re- 
sponse kernel a lJ '(t — t'). The Kubo expression for this 
response kernel is a ,J '(r) = 0(r) K l i(r). Its Fourier 
transform is the generalized susceptibility x 1 HlS). From 
here we can derive the expression (F k ) = — J^j G fc -' Xj, 
where G J is the generalized conductance matrix: 



G y = lim 



K ij (r)TdT (3) 



Following Berry and Robbins ^(J we split the conduc- 
tance matrix into symmetric and anti-symmetric parts. 
Namely G lJ = rfi + B 4 - 7 . The antisymmetric part B can 
be regarded as a vector B = (B 23 , B 31 , B 12 ), and the 
expression for the current can be written in an abstract 
way as (F) = rj ■ x — B A x. 
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The rate of dissipation, which is defined as the rate 
in which energy is absorbed into the system, is given by 
W = — (F) -x = Ylkj V^iiij ■ Only the symmetric part of 
G 13 is responsible for dissipation of energy. The adiabatic 
regime is defined by the condition |x| <C A 2 /ha, where 
A is the typical level spacing, and a is the root mean 
square value of the matrix element (dH./dx) nm between 
neighboring levels. In the adiabatic regime rf 3 vanishes 
because of the discreteness of the energy spectrum [Tn| . 
But outside of the adiabatic regime the levels acquires an 
effective width T/A = ((ha/A 2 )V) 2/3 > 1 and therefore 
the smoothed version of K lJ (ui) should be considered. 
Consequently one can obtain the fluctuation-dissipation 
(FD) relation: rf 3 ~ C %3 (uj = 0). The formulation of the 
exact FD relation depends on the assumptions regarding 
the occupation f(E n ) of the energy levels. See [ia Il4| . 
Commonly one assumes a zero temperature Fermi occu- 
pation, but this is not essential for the following anal- 
ysis. In order to derive the above expression for T we 
have used the result of 0] (Sec. 17) for the "core width" 
at the breaktime t = t prt of perturbation theory. Note 
that in the semiclassical limit (small Ti) the adiabaticity 
condition always breaks down. 

The antisymmetric part B of G 13 does not have to 
vanish in the adiabatic limit. It can be obtained from the 
adiabatic equation by looking for a first-order stationary- 
like solution Q) 01 Ell' but we prefer to regard it as a term 
in the (full) Kubo expression Eq.@. In it has been 

demonstrated that it can be written as 



B' 



-2hJ2 f{E n ) Im 
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E„ 



(4) 



Note that the "vertical" component of B vanishes in the 
"horizontal" £3=0 plane due to time reversal symmetry. 

Disregarding a possible persistent current contribution 
(that does not exist in the case of a planar $=0 cycle), 
the expression for the pumped charge is: 



Q = (t Idt = - 



T] • dx 



BAdx 



(5) 



fc=3 



If we neglect the first term, which is associated with the 
dissipation effect, and average the second ("adiabatic") 
term over the flux, then we get 



01. 



1 

2nh 



B • dx A dx — integer (6) 



The integration should be taken over a cylinder of ver- 
tical height 27tS,, and whose basis is determined by 
the projection of the pumping cycle onto the (x%,X2) 
plane. The last equality is argued as follows: The flux 
(1/fi) Jf B • dx Adx through a surface that is enclosed by 
a cycle is the Berry phase 0. The result should be in- 
dependent of the surface. Therefore the flux through a 



closed surface should equal 2ttx integer. Integrating over 
a cylinder, as in Eq.©, is effectively like integrating over 
a closed surface (because of the 2tt periodicity in the ver- 
tical direction). This means that the flux averaged Q of 
Eq.© has to be an integer. 

Before we discuss the quantum mechanical pumping, 
it is instructive to bring two simple examples for classical 
pumping. In the following we consider one particle (r) in 
a two dimensional ring as in Fig. la. 

The first example is for classical dissipative pumping. 
The conductance G = G 33 can be calculated for this sys- 
tem 01 leading to a mesoscopic variation of the Drudc 
formula. The current is / = — G x $, where — $ is the 
electro-motive-force. Consider now the following pump- 
ing cycle: Change the flux from $1 to $2, hence pumping 
charge Q = — G(l) x (<5> 2 — $1). Change the conductance 
from G(l) to G(2) by modifying the shape of the ring. 
Change the flux from $2 back to $1, hence pumping 
charge Q(2) = —G{2) x ($1 - $ 2 )- Consequently the net 
pumping is Q = (G(2) - G(l)) x ($ a - $ : ). 

The second example is for classical adiabatic pump- 
ing. The idea is to trap the particle inside the ring by 
a potential well C/ trap (ri— x\{t), r 2 — X2(t)). Then make a 
translation of the trap along a circle of radius R, namely 
x(t) = (R cos(r^), R sin(rit), $=const). It is a-priori 
clear that in this example the pumped charge per cy- 
cle is Q — 1, irrespective of <&. Therefore the B field 
must be 



B 



(xi,x 2l 0) 
2-k(x\ 



(7) 



This can be verified by calculation via Eq.Q). The sin- 
gularity along the x$ axis is not of quantum mechanical 
origin: It is not due to degeneracies, but rath er due t o 
the diverging current operator (dH/dx 3 oc l/^Jx\ + x' 2 ,)- 

We turn now to the quantum mechanical case. Con- 
sider an adiabatic cycle that involves a particular en- 
ergy level n. This level is assumed to have a degeneracy 
point at (xf\ x^\ $^ '). It follows that in fact there 
is a vertical "chain" of degeneracy points that are lo- 
cated at (x^\ x^, +2nh x integer). These degen- 
eracy points are important for the geometrical under- 
standing of the B field, as implied by Eq.QJ. Every de- 
generacy point is like a monopole charge. The total flux 
that emerges from each monopole must be 27r?ix integer 
for a reason that was explained after Eq.JBJ). Thus the 
monopoles are quantized in units of U/2. 

The B field which is created (so to say) by a vertical 
chain of monopoles may have a different "near field" and 
"far field" behavior, which we discuss below. (Later we 
further explain that " near field" means regions in x space, 
in the vicinity of degeneracy points, where » 1, while 
"far field" means regions where gx *C 1). The far field 
regions exist if the chains are well isolated. The far field 
region of a given chain is obtained by regarding the chain 
as a smooth line. This leads qualitatively to the same 
field as in Eq.Q). Consequently, for a "large radius" 
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pumping cycle in the $ = plane, we get \Q\ w 1. In the 
following we are interested in the deviation from "exact" 
quantization: If cf>^ — Owe expect to have \Q\ > 1, while 
if <^>(°) = 7r we expect \Q\ < 1. Only for the </> averaged 
Q of Eq. © we get exact quantization. 

The deviation from \ Q\ rs 1 is extremely large if we con- 
sider a tight pumping cycle around a <^ ' = degeneracy. 
After linear transformation of the shape parameters, the 
energy splitting A = E n — E m of the energy level n from 
its neighboring (nearly degenerated) level m can be writ- 
ten as A = ((x^x^) 2 + {x 2 -xf ] ) 2 + c 2 ^-^ )) 2 ) 1 / 2 
where c is a constant. The monopole field is accordingly 



, (0) I'D \V)\ 

g = ,£ (Xl-X{ ',3:2-4 ,X 3 ~x\ ) 

2 ((xi-x^) 2 + (x 2 -4 0) ) 2 + (f ) 2 (x 3 -4 0) ) 2 ) 3/2 



,(°) 



(8) 



where the prefactor is determined by the requirement 
of having a single (h/2) monopole charge. Assuming a 
pumping cycle of radius R in the $ = plane we get 
from the second term of Eq.JSJ) that the pumped charge 
is Q = TKy/gr, where g T = (<9 2 A/<90 2 )/A = c 2 /R 2 
is a practical definition for the Thouless conductance in 
this context. It is used here simply as a measure for the 
sensitivity of an energy level to the magnetic flux <3>. 

What we want to do in the following is to " interpolate" 
between the "near field" result, which is Q = 0(y/gr), 
and the "far field" result, which is Q = 0(1). For this 
purpose it is convenient to consider a particular model 
that can be solved exactly. We consider a ring with two 
barriers. The model is illustrated in Fig. 2. A version of 
this model, where the two barriers are modeled as "delta 
functions" , has been analyzed in in case of open ge- 
ometry. Below we are going to analyze a different version 
of the two barrier model, that allows an exact solution 
for closed geometry. 

We can classify the eigenstates of the closed ring into 
two categories: wire states, and dot states (Fig. 2a). The 
latter are those states that are localized in the "dot re- 
gion" in the limit of infinitely high barriers. In case of 
zero temperature Fermi occupation we define Ef as the 
energy of the last occupied wire level in the limit of in- 
finitely high barriers. The two "shape" parameters are 
the the bias x±, and the dot potential x 2 . The bias deter- 
mines whether the dot tends to exchange particles via the 
left or via the right barrier. The dot potential is loosely 
defined as the energy of the dot level (Fig. 2a). A model 
specific definition of these parameters in the context of 
the 3-site lattice Hamiltonian will be given later. 

The pumping cycle is assumed to be in the $ = plane, 
so there is no issue of "conservative" persistent current 
contribution. We start with a positive bias (xi > 0) 
and lower the dot potential from a large x 2 > Ep value 
to a small x 2 < Ef value. As a result, one electron is 
transfered via the left barrier into the dot region. Then 
we invert the bias (x\ < 0) and raise back x 2 . As a 
result the electron is transfered back into the wire via the 
right barrier. A closer look at the above scenario (Fig. 2b) 
reveals the following: As we lower the dot potential across 



a wire level, an electron is adiabatically transfered once 
from left to right and then from right to left. As long as 
the bias is positive (x\ > 0) the net charge being pumped 
is very small (|Q| <C 1). Only the lowest wire level that 
participate in the pumping cycle carries Q = 0(1) net 
charge: It takes an electron from the left side, and after 
the bias reversal it emits it into the right side. Thus 
the pumping process in this model can be regarded as 
a particular example 5] of an adiabatic transfer scheme: 
The electrons are adiabatically transfered from state to 
state, one by one, as in "musical chair game". 

For a single occupied level the net Q is the sum of 
charge transfer events that take place in few avoided 
crossings. For many particle occupation the total Q is 
the sum over the net Qs which are carried by individual 
levels. For a dense zero temperature Fermi occupation 
the summation over all the net Qs is a telescopic sum, 
leaving non-canceling contributions only from the first 
and the last adiabatic crossings. The latter involve the 
last occupied level at the Fermi energy. 

In order to get a quantitative estimate for the Q in a 
given avoided crossing, we consider the simplest version 
of the "two barrier model" that still contains all the es- 
sential ingredients: This is a three site lattice system. 
The middle site supports a single "dot state", while the 
two other sites support two "wire states". The Hamilto- 
nian is 



H 





Cl 



-i<f> 




(9) 



The three parameters are the bias x\ = c\ — c 2 , the 
dot potential x 2 — u, and the flux X3 — <I> = Ti<j). 
For presentation purpose we assume that < Ci, C2 <C 1. 
The eigenstates are E n . Disregarding the coupling be- 
tween the "wires" and the "dot" we have two wire states 
with E = ±1, and a dot state with E = u. Taking 
into account the wire-dot coupling we find that there 
are two vertical chains of degeneracies. The u f=a — 1 
chain is (0, — 1+c 2 , 2ith~ x integer) and the u rs 1 chain is 
(0, +l+c 2 , 7T + 2irh X integer). 

The eigenvalues E n are the solutions of a cubic equa- 
tion. Rather than writing the (lengthy) analytical ex- 
pressions for them we give a numerical example for their 
dependence on u in the inset of Fig. 3. The eigenstates 
are 



\n(x)) 




c 2 e 



ClC 




(10) 



where S is the normalization. Note that for E = ±1 we 
have S* = 2(c\ ± C2) 2 , while for E = we have S rs 1. 
After some algebra we find that the first component of 
the B field in the <& = plane is 



B 1 = -21m ( ^-n(x) 



Tr-n(x) 



" ( l C2 's 2 du 
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Which is illustrated in Fig. 3. For a pumping cycle around 
the u « =Fl vertical "chain" the main contribution to Q 
comes from crossing the u f=s =pl line. Hence we get 



q = ±n^£ = ±yr±2^ (ii) 

ci T C2 

where the Thouless conductance in this context is de- 
fined as gr = 2ciC2/(ci =F C 2) 2 - In both cases we have 
approximate quantization Q = ±1 + O(gx) for gr < 1, 
while for a tight cycle either Q — > oo or Q — > depend- 
ing on which line of degeneracies is being encircled. If 
the pumping cycle encircles both "chains" then we get 
Q = Ac\C2/{c\ — cf). In the latter case Q — 0{qt) for 
<7t 1) with no indication for quantization. 

For a pumping in a dot-wire system (see illustration 
in Fig. lb), in the limit of a very long wire (many sites) 
we express the Kubo formula for the conductance matrix 
using the S matrix of the dot region. The derivation 
assumes "quantum chaos" , and leads to 

G V = ^trace ( P^sA (12) 
2ni \ dXj J y ' 

This is easily identified as the Biittiker-Pretre-Thomas 



formula [llj , which has been derived for quantum pump- 
ing in open systems (Fig.le). In particular we get 
G 33 = (l/(27rft))trace(PSO-P)5f), which is just the 
Landauer formula [jjj [l(| • 

In summary we have shown how the Kubo formalism 
can be used in order to derive both classical and quantum 
mechanical results for the pumped charge Q in a closed 
system. In this formulation the distinction between dissi- 
pative and non-dissipative contributions is manifest. The 
dissipative contribution to the pumping can be neglected 
in the adiabatic regime. However, if the adiabaticity con- 
dition is violated it does not mean automatically that we 
have a dissipative effect. Classical pumping by transla- 
tion is an obvious example. For the derivation of the 
dissipative part of the Kubo formula it is essential to 
realize that in generic circumstances (unlike the case of 
translations) the adiabatic equation does not possess a 
stationary solution. 
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FIG1. Illustration of a ring system (a). The shape of the ring is controlled by some parameters x\ and X2- The flux through the ring 
is £3 = <3>. A system with equivalent topology, and abstraction of the model are presented in (b) and (c). The "dot" can be represented 
by an S matrix that depends on x\ and X2- In (d) also the flux xz is regarded as a parameter of the dot. If we "cut" the wire in (d) we 
get the open two lead geometry of (e). 




position time 
FIG2. Schematic illustration of quantum pumping in a closed wire-dot system. The net charge via the third level (thick solid line on 
the right) is vanishingly small: As the dot potential is lowered an electron is taken from the left side (first avoided crossing), and then 
emitted back to the left side (second avoided crossing). Assuming that the bias is inverted before the dot potential is raised back, only 
the second level carry a net charge Q = 0(1). 




FIG3. The first component of the B field for a particle in the middle level of the 3 site lattice model. It is plotted as a function of the 
dot potential X2 = u. The other parameters are = 0, and ci = 0.1, while C2 = 0.04 for the thick line and C2 = 0.02 for the thin line. In 
the limit C2 — > 0, all the charge that is transfered from the left side into the dot during the first avoided crossing, is emitted back into the 
left side during the second avoided crossing. Inset: The eigenenergies E n (x) for the C2 = 0.04 calculation. 



